Benzene formation in the inner regions of protostellar disks 
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ABSTRACT 

Benzene (c-CgHg) formation in the inner 3AU of a protostellar disk can be 
efficient, resulting in high abundances of benzene in the midplane region. The 
formation mechanism is different to that found in interstellar clouds and in pro- 
toplanetary nebulae, and proceeds mainly through the reaction between allene 
(C3H4) and its ion. This has implications for PAH formation, in that some frac- 
tion of PAHs seen in the solar system could be native rather than inherited from 
the interstellar medium. 

Subject headings: astrochemistry — planetary systems: protoplanetary disks 



1. Introduction 

Benzene is the smallest aromatic molecule, and is presumed to be a basic building block 
in the formation of polycyclic aromatic hydrocarb ons (PAHs). PAH s have been observed in 
protostellar disks around both low-mass T Tauri flGeers et al.ll2006l ) and intermediate- mass 
Herbig Ae/Be (HAeBe) stars (lAcke &: van den Anckerll2004l ). The number of detections in 
T Tauri-type disks is low, although this is to be expected since the incident UV radiation 
field is several orders of magnitude smalle r than in HAeBe st ars. Models indicate that the 
PAH emission comes from surface layers (IHabart et al.ll2004j ). leading to speculation as to 
whether PAHs form in these high UV environments or elsewhere. 



Both PAHs and benzene are also seen in protoplanetary nebu lae (PPNe) - PAHs in the 
circumstellar tori and in the bipolar outfl ows (IMatsuura et al.ll2004j ). and benzene in the torus 
around CRL618 (jCernicharo et al.ll200ll ). However, the processes that form benzene depend 
on the environment, with different p athways domina ting in the interstellar medium (ISM), 
in PPNe, in planetary atmospheres (jLebonnoisI 120051 ) and, as we shall show, in protostellar 
disks. In this Letter we consider the formation of benzene in a protostellar disk, and the 
implications this has for PAH formation. 
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2. Benzene formation in different environments 

The gas-phase synthesis of benzene proceeds via different reaction pathways in different 
environments, depending on physical and chemical conditions. One pervading condition for 
high abundances of benzene is that of high density (>10 9 cm~ 3 ). Several reaction schemes 
have been suggested for benzene formation: 



The interst ellar medium A re action scheme for the gas-phase formation of benzene was 
proposed by iMcEwan et al.l (119991 ) following the measurement and calculation of new rates 
for radiative association reactions. The formation pathway for benzene they suggested is 
thus: 



C 4 H+ 



H 

C4H3" + C2H2 



C 4 H+ + hu 



C 6 H 5 + H 2 



(1) 
(2) 
(3) 

c-C 6 H+ + e - — » c-C 6 H 6 + H (4) 

and variations of the above (for instance, C2H3 may replace C2H2 in reaction [2]). This 
produces moderate amounts (~10 -9 , with respect to H2) of benzene in a dense (10 4 cm -3 ) 
interstellar cloud model. 



c-C 6 H+ + hu 
c-C 6 H| + hu 
c-C 6 H 6 + H 



Protoplanetary nebulae IWoods et al.l (120021 . 120031 1 applied the above reaction scheme 
to the dense torus (n ~10 9 cm -3 ) around the protoplanetary nebula CRL618. CRL618 
i s rich in hydrocarbons , and benzene has been detected in infra-red spectra of this object 
(iCernicharo et al.ll200ll ). In the highly ionising environment around CRL618 C4HJI" is formed 
more efficiently via reactions of various ions with acetylene: 



HCO+ + C 2 H 2 
C2lTj~ + C2H2 
0213^ + C2H2 



C 4 H+ 



CO 
H 2 



C4H3 + H, 



(5) 
(6) 
(7) 



followed by reactions GHH The models of IWoods et al.l (120021 . 120031 1 produced a fractional 
abundance of benzene of ~10 -6 , in good agreement with observations. 



3. Model 



The protostellar disk model w e use will be described more fully in a forthcoming pub- 
lication, IWoods fc Willacy! (120061 ). We use a density and dust temperature profile based 
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on that of iD'Alessio et al.l (120011 ) (Fig. [I]). This is a flared accretion disk model, with a 
maximum dust grain size of 0.25 /im, similar to the size of interstellar dust. The disk has 
a mass accretion rate of M=10~ 8 M yr _1 and a surface density of £=100 gem -2 at 1AU, 
whilst the central star has the following properties: M*=0.7M Q , T*=4 000K and i?*=2R Q . 
Using these profiles, we calc ulate UV photon fluxes throu ghout the disk with the ray-tracing 
component of the model of lYorke fc Bodenheimerl (119991 ). This calculates the UV field due 
to the central star and the interstellar radiation field (ISRF) and also includes radiation 
scattering effects. At 10 AU we assume that the UV fiel d due to the c entral star is 50,000 
times the ISRF, in accordance with the observations of iBergin et al.l (120031 ). Despite the 
differing spectral shapes of the ISRF an d a typical T Tauri stellar field, with strong emission 
lines dominating the T Tauri spectrum (IBergin et al.ll2003l ). the formation of benzene in the 
very dense and well-shielded region we consider is not critically dependent on our choice of 
stellar UV field. 

Given gas densities, dust temperatures and UV photon fluxes at each point on our grid3, 
we are then able to solve the h eating and cooling balance of th e gas in a similar manner t o 
Kamp fc van Zadelhoflj (boOlh . iKamp fe Dullemondl tooj ) and lGorti fc Hollenbachl tooi ). 
The gas temperature can be underestimated by up to two orders of magnitude if it is assumed 
to be equal to the dust temperature throughout the disk, although this is less critical for 
the optically thick midplane region where benzene forms. Figure [1] shows the range of gas 
temperatures within the inner 3AU of the disk (the benzene formation region). 



We use a subset of the extensive UMIST Rate99 gas-phase chemical network (ILe Teuff et al. 



20001 ) . augmenting it with gas-grain interactions (freezeout, thermal desorption) and grain 
surface reactions so that our reaction set comprises approximately 2400 reactions amongst 
200 species. We include species with more carbon atoms than benzene to ensure that chain- 
lengthening does not result in a spurious build-up of benzene. We have selected reactions 
which are valid in the region 10-500 K; this upper limit is exceeded in the upper regions 
of the disk where temperatures can reach ~10,000K, although at such high temperatures 
molecules are destroyed. Hence we concentrate on the portions of the disk closer to the 
midplane, less than 2 scaleheights. 

The model follows the passage of a parcel of gas as it flows inwards from the outer regions 
of the disk (>35AU), and as such it traces a period in the evolutionary history of the disk. 
At the outer edge of the disk the parcel has the composition of a million- year old molecular 
cloud. The parcel then accretes onto the central star over a period of approximately 200 000 



1 We use a 35 AUxl6 AU (radial x vertical) grid with spacings of 0.5 AUx0.02 AU to model the very inner 
region of the disk. 
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years. As discussed by lWillacy et al.l (119981 ). the choice of initial abundances has relevatively 
little effect on the chemistry of most species. 



4. Results 

We find that benzene formation in a protostellar disk is efficient within 3AU of the 
central star, and fractional abundances of ~10~ 6 are achieved inside a radius of 2 AU (Fig. [2]). 
The most efficient formation mechanism for benzene in this region is due to the reaction 
between C3H4 and its ion, C3H4 , fo r which the rate has been measured to within ±25% in the 
laboratory (A;=7.48xl0~ 10 cm 3 s _1 ; Anicich et al. 1984 1^1. This is followed by neutralisation 
of the resulting c-CqB.^ ion in collisions with dust grains: 

C3H4 + C3H+ — > c-C 6 H+ + H (8) 



-C 6 Hf + grain — > #-C 6 H 6 + 



-H, 



(9) 



where g- signifies a species adsorbed onto a grain surface. 



In this reaction scheme, the adsorption energy of benzene becomes an important factor 
in determining the amount of gaseous benzene. Estimat es of this vary wildly: comparison 
of similar species in Table 4 of Hasegawa fc Herbstl j l993 ) leads to a value of 4 730 K, which 



agrees well with the experim ent of 



Arnett et al 



Jl988h using a graphite surface (4 752K). 



If one takes the approach of iGarrod fe Herbstl (120061 ). and adds the adsorption energy of 
acetylene to four times the adsorption energy of carbon and four times the adsorption energy 
of hydrogen, one arrives at the value of 7587K. T h is figu re lies in the upper range of an 
experiment on a graphite surface by lLozovik et al.l (119951 . 4 840-7 840 K). To quantify the 
effects of these different adsorption energies, we ran one model with a low binding energy of 
4 750K, and one model with a higher binding energy of 7580K. Results show a significant 
decrease in the extent of the gas-phase benzene distribution with the higher binding energy, 
as to be expected, with little change in the peak fractional abundance. 

The inner, midplane regions of our disk model are very dense (reaching ~10 14 cm -3 ), 
possibly dense enough for three-body (termolecular) reactions to have an effect on the chem- 
istry. In order to investigate the effect on benzene production in particular, we ran a third 
model including a representative three-body reaction: 



C3H3 + C3H3 



M 



C 6 H 6 + M, 



(10) 



2 The UMIST reaction network does not make a distinction in this case between isotopomers of C3H4; 
however, the experimental data for reaction ((8]) was obtained using allene, rather than methylacetylene 
(CH3CCH) or cyclopropene (C-C3H4) 
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where two propargyl radicals collide with M, a third molecule. This has been shown to be a 
ve ry efficient forma t ion m echanism for benzene when densities are high; see the discussion 
of ICherchneff et al.l (Il992l ). We have adopted a rate coefficient from chemical models of the 
Saturnian atmosphere: 

Kter = l — T CIIl's , (11) 

koon + k 

(IMoses et al.l 120001 ) where , the high pressure limit to the reaction rate, was taken t o 
be 1.66xl0~ 13 cm 6 s _1 and k , the low pressure limit, lxlO -27 cm 3 s _1 (jWong et al.ll2000l ). 
Inclusion of this reaction in the model did not result in a significant change in the abun- 
dance of gas-phase benzene, despite the presence of large fractional ab undances o f C3H3 
in dense regions. Even adopting the faster rates for koo and k from iLebonnoisI (120051 ) 
(1.2xl0 _10 cm 3 s _1 and 6.0xlO~ 28 e 1680//T cm 6 s -1 , respectively) had little overall effect. 



5. Discussion 



The region in which benzene is abundant in our model (Fig. [2]) is also rich in other 
organic molecules: high densities mean fast collision timescales with other molecules and 
with dust grains, and also protection from the strong stellar UV field in this region. These 
prime conditions for the growth of organic molecules may allow molecules to grow beyond the 
single arom atic ring of benzen e into PAHs, either by substitution of hydrogen by acetylene 
in the ring ( jWong et al.l l2000l : IMoses et al.l 120001 ) . or by the dimerisation of benzene (Ai) 
with successive phenyl rings to produ ce naphthalene (A 2 ), phenanthrene (A3), pyrene (A 4 ) 
and successively large PAH molecules (I Cher chnefB 1 1 9 96l ) . Benzene itself is unlikely to persist 
through to the formation of planets: calculations show that benzene has a timescale for 



UV destruction in the diffuse solar system of only hundreds of years (lAUain et al.l 11996 



Ruiterkamp et al.l 120051 ) . meaning that the benzene seen on Jupiter, Saturn and Titan is 
unlikely to be primordial. 

In order for PAHs to form, the constituents must be present in a region warm enough 
for ring closu re to occur. PAHs only seem to form efficiently at temperatures in the region of 
900-1 100 K ( IFrenklach fc Feigelsonlll989l ) and not greater. Production may be possible, too, 
at lower temperatures (700-900 K) given a high pressure (jHelling et al.l 119961 ) . The region 
of benzene production is at a temperature of a few hundred degrees, and sits just below a 
region of sufficiently high temperature for PAH formation. Invoking vertical mixing here 
could raise benzene and other organics into the higher temperature region, as well as push 
newly-formed PAHs into the surface layers of th e disk, where they ar e observed. Mixing 
timescales at a few AU are fast enough (~625yr; lllgner fc Nelsonll2006l ) to dredge material 
from the midplane to surface layers within the lifetime of a T Tauri disk (~10 6 yr). 
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For PAHs to subsist in the solar system, the rate of growth of the PAH must be greater 
than the rate of photodissociation. However, once a PAH reaches 30-40 carbon atoms in 
size, it is effectively safe from photodissociation since at this size the infrared radiative 
rate of the PAH dominates the photodissociation rate. In the ISM, a benzene molecule is 
destroyed (by removal of an acetylene molecule) at a rate of 1.5 xlO -10 s -1 , giving a lifetime 
of around 200 years. The growth time (acc retion of an acetylene molecule) is given by 



r gr = (n X C2i i 2 kc 2 H 2 ) 1 s ( jAllain et al.lll996l ). Inserting typical values for the abundance 
of acetylene in the inner part of the disk, X C2 h 2 =10~ 7 , and the reaction rate for accretion 
of acetylene, &c 2 h 2 = 10 -11 s 5 gives a necessary density of 1.5xl0 8 cm~ 3 for PAH growth 
from benzene. This tallies with the low abundances of benzene in (low density) interstellar 
cloud models (see Sect. [2]) and the lack of a detection of benzene in the ISM. However, these 
conditions are met within the inner regions of a protostellar disk, as has been shown in our 
chemical model (which also includes destruction of benzene by reactions with ions). In fact, 
the strength of the UV field drops below the interstellar level for much of the midplane region, 
and will remain at that level for much of the transition towards a debris disk. Assuming a 
PAH of 40 carbon atoms forms solely through the addition of acetylene to a benzene ring, a 
PAH could be produced from a benzene molecule on a timescale of approximately fifty years 
at a density of 10 10 cm -3 . Even though the actual route to PAH formation is likely to be 
different, this estimate does give a representative timescale and show that PAH production 
from benzene could possibly be prolific in this type of environment. We will consider the 
formation of larger aromatics in a future paper. 

Having ascertained that benzene can survive for an amount of time long enough for 
PAHs to form from it, the question arises, do PAHs form in similar regions to which benzene 
forms? The answer to this is unclear, since the suppor ting observational evidence is sparse. 



Of the three detections of PAHs in T Tauri-type disks (IGeers et al.ll2006l ). all three disks are 
thought to have inner dust holes (from SED models). Thus it may just be a selection effect 
that PAH emission is only detected from the upper layers of disks because of the presence 
of these dust holes: otherwise the continuum emission swamps any possible PAH emission 
feature originating from deeper in the disk. Further observations are necessary to confirm 
the correlation between the presence of benzene and PAHs. 



This research was conducted at the Jet Propulsion Laboratory, California Institute of 
Technology under contract with the National Aeronautics and Space Administration. Sup- 
port comes from an appointment to the NASA Postdoctoral Program at JPL, administered 
by ORAU through a contract with NASA. 
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Fig. 1. — Density and gas temperature profiles for the inner 3AU of the disk model. The 
left panel shows gas density, whereas the right panel shows the gas temperature calculated 
by the heating/cooling balance. See the electronic edition of the Journal for a colour version 
of this figure. 
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Fig. 2. — The inner 3 AU of the disk, showing the distribution of (gaseous) benzene in terms 
of fractional abundance. The black solid lines are optical depth (r) contours. The left panel 
shows a model with a low binding energy for benzene, 4 750 K. The right panel shows the 
effect of a higher binding energy, 7580K. See the electronic edition of the Journal for a 
colour version of this figure. 



